Observation of the decay B° — > J ftp r\ 
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We report the first observation of B° — ► J/ip 77 decay. These results are obtained from a data 
sample that contains 449 million BB pairs accumulated at the T(4S) resonance with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. We observe a signal with a significance of 
8.1ct and obtain a branching fraction of (9.5 ± 1.7 (stat) ± 0.8 (syst)) x 10 -6 . 
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At the quark level, the decay B° — ► J/tpr) proceeds 
primarily via a6-> ccd transition. As is apparent from 
its Feynman diagram (Fig. [I}, this is a Cabibbo- and 
color-suppressed decay, similar to B° — > J/ipir° for which 
the branching fraction and CP violation parameters have 
been measured by the BaBar and Belle collaborations [l], 
& Si EL H @] ■ The best existing limit for the branching 
fraction for B° — > J/iprj comes from BaBar pj and is 
based on 56 million BB pairs. 

If factorization and flavor-SU(3) symmetric coefficients 
for the dd content are assumed, the branching fraction for 
B° — > J/ipr/ decay is expected to be comparable to that 
for B° — * J/i[m°- If the measured branching fraction is 
significantly different from this expectation, it would im- 
ply the influence of large final-state interactions or non- 
standard contributions. In the latter case, precise mea- 
surements of CP violation parameters might also reveal 
unexpected phenomena. 

In this letter, we report the first observation of the 
decay B° — > J/t/jt] The measurement is based on a 
data sample that contains 449 million BB pairs, collected 
with the Belle detector at the KEKB asymmetric-energy 
e+e- (3.5 on 8 GeV) collider [| operating at the T(45) 
resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a barrel- like ar- 
rangement of time-of- flight scintillation counters (TOF), 
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FIG. 1: Feynman diagram for the leading contribution to the 
decay B° — > J/ip n. 

and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL). All these detectors are located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect K® mesons and to iden- 
tify muons (KLM). The detector is described in detail 
elsewhere [lOj. 

The data used in this analysis was collected with two 
detector configurations. A 2.0 cm radius beampipe and 
a 3-layer silicon vertex detector are used for the first 
sample of 152 million BB pairs, while a 1.5 cm radius 
beampipe, a 4-layer silicon detector and a small-cell in- 
ner drift chamber are used to record the remaining 297 
million BB pairs [ll|. A GEANT-based Monte Carlo 
(MC) simulation is used to model the response of the 
detector and determine efficiencies flil ]. 

Charged tracks are selected based on the impact pa- 
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ramcters relative to the interaction point: dr for the ra- 
dial direction and dz for the direction along the positron 
beam. Our requirements are dr < 1 cm and \dz\ < 5 
cm. From among the selected charged tracks, e + and e~ 
candidates are identified by combining information from 
the ECL, the CDC {dE/dx) and the ACC. Identification 
of fi + and candidates is based on track penetration 
depth and the hit pattern in the KLM system. Charged 
pions are identified using the combined information from 
the CDC (dE/dx), the TOF and the ACC. 

Photon candidates are selected from calorimeter show- 
ers not associated with charged tracks. An energy de- 
position with a photon-like shape and with energy of at 
least 50 MeV in the barrel region or 100 MeV in the end- 
cap region is considered a photon candidate. A pair of 
photons with an invariant mass 117.8 MeV/c 2 < M 77 < 
150.2 MeV/c 2 is considered as a ir° candidate. This in- 
variant mass region corresponds to a ±3cr interval around 
the 7T° mass, where a is the mass resolution. 

We reconstruct J/tp mesons in the £ + £~ decay channel 
(i = e or ji) and include up to one bremsstrahhmg photon 
that is within 50 mrad of each of the e + and e~ tracks (de- 
noted as e + e~(j)). The invariant mass is required to be 
within -0.15 GcV/c 2 < M ee(7) - m J/v , < 0.036 GcV/c 2 
and -0.06 GcV/c 2 < M^-mj/y, < 0.036 GeV/c 2 , where 
mj/^ denotes the nominal mass, M ee ( 7 ) and M MA1 are the 
reconstructed invariant mass from e + e~(7) and ji + ji~ , 
respectively. Asymmetric intervals are used to include 
part of the radiative tails. 

Candidate rj mesons are reconstructed in the 77 and 
7r + 7T~7r° final states. We require the invariant mass to be 
in the range 500 MeV/c 2 < M 77 < 575 MeV/c 2 (±3cr) 
and 535 MeV/c 2 < M^-^o < 560 MeV/c 2 (±3cr). In 
the 77 final state, we remove rj candidates if either of the 
daughter photons form a 7r° candidate when combined 
with any other photon in the event. Asymmetric decays 
are removed with the requirement \E^\ — E^/ (E^i + 
-E72) < 0.8, where E 1 \ and S 7 2 are the energies of the 
two photons that form the 77 candidate. 

We combine the J/tp and r\ to form B mesons. Sig- 
nal candidates are identified by two kinematic variables 
defined in the T(45) center-of-mass (CM) frame: the 
beam-energy constrained mass, M^ c = ^/i? 2 eam — p* B 2 , 
and the energy difference, AE = Eg — -Ebeamj where p* B 
and E* B are the momentum and energy of the B can- 
didate and -Ebcam is the run-dependent beam energy. 
To improve the AE resolution, the masses of the se- 
lected 7T°, rj and J/tp candidates are constrained to their 
nominal masses using mass-constrained kinematic fits. 
In addition, vertex-constrained fits are applied to the 
r) — > 7T + 7r~7r° and J/tp — > candidates. We retain 

events with M hc > 5.2 GeV/c 2 and \AE\ < 0.2 GeV, 
and define a signal region to be 5.27GeV/c 2 < M^ c < 
5.29GeV/c 2 , -0.10 GeV < AE < 0.05 GeV (rj -► 77) or 
\AE\ < 0.05 GeV (77 — > Tr + n~ir a ). In events with more 
than one B candidate, usually due to multiple rj candi- 



dates, the B candidate with the minimum x 2 value from 
the mass- and vertex-constrained fit is chosen. 

To suppress the combinatorial background dominated 
by the two-jet-like e + e~ — * qq (q = u, d, s) contin- 
uum process, we remove events with the ratio of sec- 
ond to zeroth Fox- Wolfram moments R2 > 0.4 [l3j |. 
This requirement is determined using a figure-of-merit 
Ns I \AZVg + Nb, where Ns is the number of expected sig- 
nal events and Nb is the number of background events. 
For Ns, we use a MC simulation with the assumption 
that B(B° -> J/tp if) is 6 x 10" 6 [3]. For N B , we use a 
sample of continuum background MC that corresponds to 
an integrated luminosity that is about three times that of 
the data sample and normalized to the number of events 
in an M hc - AE sideband (5.2 GeV/c 2 < M hc < 5.26 
GeV/c 2 and 0.1 GeV < AE < 0.2 GeV). This require- 
ment on i?2 eliminates 88% of the continuum background 
and retains 97% of the signal. 

After the above continuum suppression, the back- 
ground is dominated by BB events with B decays to 
J/tp. We use a MC sample corresponding to 3.86 x 
10 10 generic BB decays that includes all known B° — > 
J/tpX processes to investigate these backgrounds. We 
find that the dominant backgrounds come from B° — > 
J/ipK L (25.1%), B ± -> J/tpK* ± (892) (19.3%), B° -> 
J/il)K*°(m2) (14.5%), B° -> J/i'Ks (14.1%), B° -> 
J/ipw (8.4%), and a few other exclusive B — * J /if) X 
decay modes. These backgrounds peak in the Mbc distri- 
butions but not in AE. Thus, combinatorial backgrounds 
from B — > J/tp X decays do not affect the signal yield 
extracted from the AE distribution. 

Signal yields and background levels are determined by 
fitting the distributions in AE for candidates in the Mb c 
signal region. The AE distribution is fitted with a func- 
tion that is the sum of a Crystal Ball line function [l5| 
and two Gaussians for signal, and a second-order polyno- 
mial for background. The probability density functions 
(PDF) that describe the signal and background shapes 
are determined from MC simulations. We use the MC 
sample for this process to determine the PDF for the 
total background. 

We use the decay B+ -> J/ipK*+(K*+ -> K+ir°) as 
a control sample to correct for differences between data 
and MC for the fitted mean and width values of the AE 
signal peak. For this selection, we require the helicity 
angle in the K* + — > K + ir° decay, i.e. the angle between 
the 7T° and the negative of the B + momenta in the K* 
rest frame, to be less than 86 degrees. This requirement 
primarily selects events with a high momentum 7r° and 
produces a control-sample AE distribution that is similar 
to that for our signal decay mode. The signal PDF is 
modified based on the results for the control sample: the 
mean value of AE is shifted by -3.81 ± 0.65 MeV and 
the width is scaled by a factor of 0.99 ± 0.04. 

There are 98 and 58 events in the Mb c signal region 
for the rj — > 77 and r) — > n + ir~ir° modes, respectively. 
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We determine the signal content by performing an un- 
binned extended maximum-likelihood fit to the candi- 
date data events. The unbinned extended maximum- 
likelihood function is 

(Ns+N B ) JL. 

£=— ^ Y[[N s P s {AEi) + N B PB(AEi)} , (1) 



where N is the total number of candidate events, N$ 
and Nb denote the number of signal and BB back- 
ground events, Ps(AEi) and Ps{AEi) denote the sig- 
nal and background AE PDFs respectively, and i is 
the event index. Separate fits to the B° — ► 7/^77(77) 
and B° — > J / iIj7j(tt + tt°) candidate samples give sig- 
nal yields of 43. f ± 8.9 and 16.6 ± 5.8 events, respec- 
tively. The results of the fits to the data are shown 
in Fig. O These signal yields correspond to branching 
fractions of (9.6 ± 1.9) x 1CT 6 and (10.0 ± 3.5) x 10~ 6 , 
respectively. We fit the two samples simultaneously, con- 
straining the results to a common branching fraction, and 
obtain 59.7 ± 10.5 events. The combined branching frac- 
tion (9.5 ± 1.7) x 10~ 6 . 

We correct for discrepancies of the signal detection ef- 
ficiency between data and MC using control samples. 
The J/ip -> and D*+ -> D°tt+ (D° -> R--K+) 

events are used as control samples to correct for the 
muon identification and pion identification, respectively. 
We find that the MC overestimates the efficiencies for 
J/ip ->■ [i+ii- and 7/ -> tt+tt-tt by 8.5% and 1.8%, re- 
spectively. 

A 4.8% systematic error associated with uncertainties 
in the signal and background PDFs is estimated by com- 
paring the fit results for the cases when the polynomial 
parameters are fixed by either MC or data in the Mb c 
sideband region, and from changes that result from vary- 
ing each parameter by one standard deviation. For the 
contribution from the corrections to the lepton identifica- 
tion efficiencies, we use a control sample of J/ip — > £ + £~, 
which indicates uncertainties of 2.7% per electron track 
and 1.2% per muon track. By weighting the relative con- 
tributions from J/ip — > e + e~ and J/ip — » fi + \i~ , we 
assign a 3.9% systematic error. The systematic error 
from the pion identification correction is 0.7% per track; 
this is determined from a study of the D* + — > D°tt + 
(D° — > K~ir + ) control sample. Since this applies only 



to the 77 



~7T° decay mode, the effective systematic 



error is 0.4%. The systematic error due to the track find- 
ing efficiencies for the Belle detector is 1.0% per charged 
lepton and 1.3% per charged pion. This contributes 2% 
per J/ip and 0.7% per 77, and the total of 2.7% (from lin- 
early adding) is included. A 4.1% systematic error due 
to 7T° detection is determined from a comparison of the 
data and MC ratios for a large sample of 77 — ► 7r + 7r~7r° 
and 77 — > 37T° decays. Since 77 — ► 77 is similar to 7r° decay, 
we also assign 4.1% as the systematic error for 77 — > 77 
reconstruction. The errors on the branching fractions for 
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FIG. 2: AE distributions for the decay B° -> J/ip-q. The 
plot on the top shows B° — > J/ipr](i] — > 77) while B° — > 
J/tprj(r} — > 7T + 7r - 7r ) is on the bottom. The curves show the 
signal (dash-dot lines) and background (dashed lines) contri- 
butions as well as the overall fit (solid lines). 



J/V -> e+e~, 77 -> 77, and 77 -> tt+tt-tt are 1.0%, 0.7%, 
and 1.8%, respectively [T3|. These errors contribute 1% 
and 1.1% systematic errors due to B(J/tp £ + £~) and 
/3(t7 — * 77 and 7T + 7r~7r ), respectively. We assign a sys- 
tematic error of 1.2% due to the uncertainty in the num- 
ber of BB pairs. 

The systematic errors are summarized in Table HI The 
total uncertainty is the quadratic sum of each term. The 
detection efficiencies and branching fractions are listed 
in Table |TTJ The statistical significance of the observed 
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signal, defined as \J — 2 ln(£o /^ma) where £ max (£o) de- 
notes the likelihood value at the maximum (with the sig- 
nal yield fixed at zero)- is 8.1a. 

TABLE I: Systematic errors for the combined branching frac- 
tion. 



Systematic errors(%) 



PDFs 


4.8 


Lepton identification 


3.9 


Charged pion identification 


0.4 


Tracking (lepton and charged pion) 


2.7 


rj — > 77, 7T° selection 


4.1 


B(J/ip ->i+r) 


1.0 


B(rj — > 77 and 7r + 7T~7r°) 


1.1 


N BS 


1.2 


Total 


8.1 


TABLE II: Detection efficiencies and branching fractions. 


Mode Eff.(%) 


BF(10~ 6 ) 


B° -> J/ipvill) 21.9 


9.6 ± 1.9 


B° -> J/^Itt+tt-tt )! 13.8 


10.0 ± 3.5 


_B° — > J/ipr] (combined) 


9.5 ± 1.7 ± 0.8 



In summary, the first observation of the decay B° — > 
J/V>f? is reported. We observe 59.7 ± 10.5 signal events 
with 8.1cr significance in a 449 million BB pair data sam- 
ple at the T(45) resonance. The measured branching 
fraction is B(B° -> J/ip rj) = (9.5 ± 1.7 ± 0.8) x 10~ 6 , 
where the first error is statistical and the second is sys- 
tematic. The measured branching fraction is consis- 
tent with the theoretical expectation and compa- 
rable to the world-average value for the J/ipir mode: 
(2.2 ± 0.4) x 10" 5 0. There is no indication of either 
large final state interactions or non-standard model con- 
tributions. 
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